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ABSTRACT 



Here to represent the propagation of waves I attempted to describe them 
separately in space and time domains. The time and space wave equations are 
obtained and investigated, and formulas for the wave propagation are 
expressed. I also tried to apply the model to the description of such physical 
phenomena as the constancy of the speed of light and propagation of ocean 
waves - tsunami. 
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1. Model Assumptions 



For generality, I will consider here 'the energy's waves'. It gives me the 
possibility to abstract from some insignificant details (from the model's point of 
view). 

For simplicity, I will describe the one-dimensional space. One can look at 
the coordinate as a vector to get a picture in multi-dimensional space. 

I denote u(x,t) to be the value of energy in point x at time t. I also 

denote <&(x,t) to be the level of energy in point x at time t (coincides with 

u(x,t) when the system is at rest). 

Difference between the value u(x,t) and the level <&(x,t) constitutes the 
energy's disturbance AU(x,t), i.e. 

AU(x,t)=U(x,t)-<t>(x,t). (1) 

When system is at rest, we can say 

AU(x,t) = (2) 
or equivalently 

U(x,t)=<t>(x,t). (3) 

I don't mean in (3) that system has to be static to reside at rest. We will 
consider a moving system in the eighth section of paper. 

In contrast to the conventional wave equations, which are described by 
partial differential equations of second order (see [3, 4, 6, 7]), I will be using the 
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following assumptions in the attempt to describe the dynamics of energy's 
disturbance separately in time and in space. 
In the time domain, 

(1 ) The second derivative of energy's value with respect to time is inversely 
proportional to energy's disturbance. 

(2) The first derivative of energy's level with respect to time is directly 
proportional to energy's disturbance. 

In the space domain, 

(3) The second derivative of energy's value with respect to direction is 
inversely proportional to energy's disturbance. 

(4) The first derivative of energy's level with respect to direction is directly 
proportional to energy's disturbance. 

To simplify the reasoning I also introduce the dual variables P(x,t), which 
I will describe in the next two sections. 

2. Model Description 

Until time t Q = I assume that system was at rest, i.e. 



U(x,t)=®(x,t)=<$ 



(4) 



for 



oo<t<0 and for all x. 



At time t = the system was excited in the point x = . 



Thus we can write 



t/(0,0)=<& +A<I> . 



(5) 



4 



Let me try to describe the dynamics of disturbance 

AU(x,t)=u(x,t)-<i>(x,t), 

t A „/ nn \ A ^ dAU (0,0) A A dAU{0,0) , 
for 0<t<+°°, 0<x<+°°, At/ (0,0)= A$ , ^-^ = o, and -^-^ = o 



dt 



dx 



In the time domain, 



^i2 t l = -^(U( Xl ,t)-^(x 1 ,t)), 



dt 



d 2 U(x l ,t)_ dP(x l ,t) 

2 A 2 



dt' 



dt 



d$>{x y ,t) _ . dP{x x ,t) 

— —A, , 

dt dt 

where ^>0, a,,a 2 ,a 3 >0 are constants, and Pfaj) is dual variable, which 
connects the second derivative of u{x x ,t) and first derivative of <&(x v t) with 
disturbance Auix^t). 

We can solve equations (6) - (8) (see [1 , 5]) to obtain the dynamics of 
disturbance AC/(x 1 ,?) in the time domain, 

d 2 d 

— - AU (xj ,t)+ A l A 3 — AU (xj ,t)+ \k 2 AU (x l , t) = . 
dt dt 



(a) If i > p- the solution of (9) is 

2 Al Aj 



(6) 
(7) 
(8) 



(9) 



At/(jc 1 ,0=Qe* 1 ' +C 2 e* 2 ', 
where Q and C 2 are constants of integration, and 



(10) 



k l 2 - Aj 



A, 



2-k 



v 2 y 



4l 
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Since both <0 and k 2 <0, the value A[/(x l5 /)^0 for ?->+oo, and 
oscillations don't take place in the time domain. 

(b) if 4. = Bl the solution of (9) is 

2 v i 

Af/(^,/ 1 )=(C 1 +C 2 / 1 )e (11) 
Again the value AC/(x 1 , * ) -> for * -> +oo , and oscillations are not here in 
the time domain. 

X 

(c) If — < J— the solution of (9) is 



^1 ^3 

AC/^ , t ) = e"^ ' (C l cos y^r + C 2 sin ) , (1 2) 



2 2 



where B = X,, 

,v 4 



Here the value AC/fo , * ) -> for ? -> +oo , and we have the situation of 

damped oscillations in the time domain. 

Now let me do the same for energy's disturbance in the space domain, 

dP{x,t x ) 



dx 



d 2 U(x,t l ) _ dPjx^) 

dx 2 ~ Ml dx ' [ ' 

Li . V Li .X 

d&(x,t l ) _ dPjxj,) 
dx dx 

where t x > , n x , fi 2 , /x 3 > are constants. 

Then dynamics of disturbance AC/(x,^) in the space domain is, 
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d 2 d 

— T-AU(x,t l )+ — AU(x,t l )+ ju l ju 2 AU(x,t l )=0. 
dx dx 



(d) if > fe. the solution of (1 6) is 

2 \Mi 



(16) 



AU(x,t l )=C l e' 1 ' +C 2 e' 2 ' . 



(17) 



2 "II 



00 

V 2 , 



0i 



where / u = ^ 

Since both ^ < and l 2 < , the value AC/(x, ^ ) -> for x -> +oo , and 



oscillations don't take place in the space domain. 



(e) If ^ = 1^ the solution of (1 6) is 

2 Mfi, 




_thth x 

AU(x,t l )=(C 1 +C 2 x)e 2 . 
Again the value AU(x,t 1 )^>0 for x 
the space domain. 



(f) If ^ < p^- the solution of (1 6) is 

2 \Mi 



Mi Mi x 

AL/(x,t l )=e 2 (C l cosyx + C 2 sin fx), 
where y = JU\- 



(18) 

and oscillations are not here in 



(19) 









V 2 J 



Here the value AU{x,t x )^0 for x ->+<*>, and we have the situation of 
damped oscillations in the space domain. 
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Earlier I said about intention to represent the dynamics of disturbance 

AU(x,t) from initial values AC/(0,0)= A<D , ^M = o,and dAU (°^ = . 

dt dx 

Let me show how the value AU(x x ,t x ) is obtained for x x > and t x > . 
If I denote c the velocity of disturbance traveling, then disturbance will 

reach the point x x at time t = x y c . 

Hence we have AU (x x , t x ) = for t x < t Q . 

For t x = t Q +At (where A t > ) we may find at first AU (0, A t) from 

equation (9), and initial values At/(o,0)= A4> and ^ a ^( Q '°) = q yh en we fj nc | 

J? 

the value AU{x x ,At) from equation (16), and initial values At/(0,A?) and 

dAUp 1 At)_ Q 

dx 

Thus the obtained value AU(x x ,A t) is the energy's disturbance in the 
point x x at the time A t for ideal 'instant' velocity of disturbance traveling that is 
equal to the value of energy's disturbance in the point x x at the time (t +At) for 
disturbance traveling with the real velocity c. 

Note we used above the boundary condition dAu ^'^ = dAU(0,0) _ Q ^ 

dx dx 

which is easy to understand through the energy 'invariant' described in the fourth 
section. 
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3. Dual Model 



Here we describe dynamics of dual variable P(x,t). 

At first in the time domain let me differentiate (6) and integrate (7). It gives 

us, 

d 2 P{x l ,t) + ^ dP ^ t ) + x i X 1 P{x l ,t)+C,=Q, (20) 



dr dt 



where C 3 = X x 



J 



. Note that the values du ^ f ^ and 

dt 



dt 

P(x v t Q ) are taken at time t when energy's disturbance reaches the point x x . 
We can simplify (20) by changing the variable, 

P 1 (x 1 ,t) = P(x 1 ,t) + ±- dU(x ^ - P ( Xl ,t ). 
A 2 dt 

It will give us, 

+ K X, + A l A 2 P 1 {x 1 ,t)=0. (21) 

dt dt 

Equation (21) can be solved similarly to (9). Hence we obtain P^,?)-* ° 

for ?^>+<>o,and therefore P(x r ,t) -» P(x l ,t Q )-— dU ^ Xl,t ^ when ? ->+<*> for all 

X 2 dt 



x x > . 



Similarly in the space domain, 

d 2 P(x,t,) dP(x,t,) , x „ rnn , 

dx dx 
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where C 4 = ju x 



r dU(x t) ,t.) 



dx 



. Note that the values ^(vO and 



J 



dx 



P{x Q ,t x ) are taken either (a) at the time t l in the excited point x Q for energy's 



disturbance traveling with ideal 'instant' velocity, or (b) at the time 



in 



the excited point x Q for energy's disturbance traveling with real velocity c. 



Therefore P{x,t x ) -» P(x Q ,t x ) - — dU ( ' X ° ,t ^ when x ->+<*> for all t x >0. 

H 2 dx 



4. Energy Invariant 



We can see that oscillations involve transformations of energy's value into 
energy's 'velocity' in the time domain and energy's 'tangent 'in the space domain. 
Therefore we can deduce the following energy's 'invariant'. 
In the time domain, 



e(x l ,t)=^(AU(x 1 ,t)) 2 + 
and in the space domain, 



~P ~dt 



kU{x x ,t) 



(23) 



e (x, t x ) = (At/ (x, t x )) 2 + — — At/ (x, ^ ) 

y dx 



1 J 



(24) 



Here I mean that if 'damping' coefficients A 3 and ju 3 were zeroes then 
both values e(x l ,t) and £(x,^) would be constants. 

Then in both time and space domains energy 'invariant' would look like, 
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e{x,t)= 



1 d 



(AU(x,t)f+ ^-~AU(x,t) 
p at 



+ 



y dx 



AU(x,t) 



(25) 



In other words, without 'damping' coefficients we would have infinite 
oscillations both in time and in space. 



5. Interference of Energy's Waves 

I will consider the situation where the system was initially at rest, and later 
it was excited twice. 

First time it was done in point x Q = at time t Q = 0, where A£/(o,0)= A<£ , 

dAU (0,0) A . dAU (0,0) A c ■ ■ ♦ ' n ♦ ♦■ 
= , and = . Second time it was in point x Q = p > at time 

dt dx 
t' Q =x>0.\ assume that system was excited second time before the energy's 

first disturbance reaches the point x' , therefore r<^/, AU(p,t)= A<£' , 

^• , ) = 0,and^- t) = 0. 

dt dx 

Let me track the value AU(x 1 ,t 1 ) for x 1 > p and t l >t . 
We have AU(x 1 ,t 1 )=0 if t x <r + ^^-. 



For 



f x x -p^ 

T + — — 

C J 



<t x < — , we will find the value AU(x 1 ,t 1 ) from initial 



values AU(p,t), dAU (P> z \ an( j dAU(p,z) s j m j| ar |y as it was done in the 

dt dx 
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second section. We only have to apply here different time interval, 

At = t l -T- — — — . 

c 

For t x > — , the energy's disturbance from both points x = and x' = p 

c 

will reach the point x 1 . At first, we may obtain the value AU^x^tJ from initial 

values AU(p,z), dAU (P> z \ anc | dAU(p,z) ^ usjn g the tjme j nterva | 

dt dx 

At l =t 1 -T-— — — ). Then we can find the value AU 2 (x l ,t 1 ) from initial values 

c 

KTT(nn \ dAU(0,0) . dAU(0,0) . . ,. ,. . , .. . x x . 

AU (0,0), and (where the time interval is At 2 =t l -—). 

dt dx c 

Finally we calculate the united value At/fa,^) by adding both partial values 
together, i.e. AU(x 1 ,t 1 )= AU l (x l ,t l )+ AU 2 (x l ,t l ). 



6. Generalization of 'Boundary' Conditions 

We have described so far the situation when system was excited in some 
point x Q with non-zero energy's disturbance AU(x ,t Q )*0 and zero disturbance 

ra{e dAU(x ,t ) = Q 
dt 

Assuming the non-zero disturbance rate in initial point x =0 and time 
t =0 doesn't bring big changes in our previous reasoning. We have to consider 



additionally the initial value d AU(x ,t ) _^ Q w h en we are solving equation (9) in 

dt 

the second section. 

Let me look now at case when the system is excited in some interval \Q,a] 

at initial time t = 0. We may assume the regular condition with AC/(0,0)= A<£ , 

dAf/(0,0) . dAt/(0,0) . , \ j T ( f\\ r\ dAU(a,0) . , dAU(a,0) . 

^_^0, ^-^*0 and AU(a,0)=0, ^-^o, and ^j_^ = o. 

dt dx dt dx 

Other points in the interval (0,a) have initial values somewhere in the middle. 

When we calculate the value Al/(x l ,t l ) for > a and t x >0, the 
interference of energy's waves has to be taken into account. 

x — a 

Hence, for t l < — we have, 

c 

AE/fo.OsO. (26) 
Interference of energy's waves appears for bigger *, . To simplify its 
representation I denote A^{x,t,p,x) the value of energy's disturbance AU(x,t), 

which spreads to the point x and time t from initial values AU(p,t), dAU (P> z \ 



and 



dt 

dAU{p,t) 

dx 



x — Cl X 

Then for — < t l < — we can write, 

c c 



a 

AU(x 1 ,t l )= \Ay(x l ,t l ,ZWZ- (27) 



For t x > — it takes place, 
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a 

AU(x 1 ,t 1 )=\A*¥(x 1 ,t 1 ,£,0)dZ- 





(28) 



7. Disturbance Traveling and Wave Propagation 

We have to distinct the concepts of velocity c of disturbance's traveling 
and velocity of wave propagation. 

The former is an external value for our model while the latter one is 
determined inside the model. 

Indeed the velocity w of wave propagation is the ratio of the period of 
oscillations in space, and the period of oscillations in time (or inverse ratio of the 
corresponding frequencies). 



T. 



space 



I 
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Thus w = 



time 



Hence 



w = — 




(29) 



and when 'damping' coefficients A 3 =0 and ^ 3 =0, we have, 




(30) 
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8. Example One: Invariance of the Speed of Light 

I want to underline again that velocity w of propagation of energy's 
disturbance is not the same as velocity c of traveling of energy itself. 

This point becomes important when we address the phenomenon of the 
speed of light. 

If me move the system associated with energy's level <i>(x,t) with non- 
zero velocity then the velocity w of wave propagation will remain the same in 
new system, and will be described by equation (29). 

Here I see the reason of negative outcome of Michelson-Morley 
experiment, which started the creation of the theory of relativity [2, 3, 4]. 

9. Example Two: Tsunami or 'Harbor Wave' 

It looks possible to apply the model to describe the phenomenon of 
'tsunami'. 

An earthquake passes its energy to ocean water. This energy is 
transformed into elevating of water column. 

Using equations (17) - (19) we can obtain the value MJ{x,t x ) of energy's 

disturbance traveling on significant distances from the earthquake's place. 

Then we can calculate the height, on which the water column of particular 
depth has to be elevated to convey the required amount of energy. 

That value gives us the height of tsunami. 
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